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ABSTRACT 

We study the dynamical stability of an additional, potentially habitable planet 
in the HD 47186 planetary system. Two planets are currently known in this sys- 
tem: a "hot Neptune" with a period of 4.08 days and a Saturn-mass planet with 
a period of 3.7 years. Here we consider the possibility that one or more un- 
detected planets exist between the two known planets and possibly within the 
habitable zone in this system. Given the relatively low masses of the known plan- 
ets, additional planets could have masses <1OM0, and hence be terrestrial-like 
and further improving potential habitability. We perform A^-body simulations 
to identify the stable zone between planets b and c and find that much of the 
inner habitable zone can harbor a 10 M0 planet. With the current radial velocity 
threshold of ~ 1 m/s, an additional planet should be detectable if it lies at the 
inner edge of the habitable zone at 0.8 AU. We also show that the stable zone 
could contain two additional planets of 10 each if their eccentricities are lower 
than ~ 0.3. 

Subject headings: stars: planetary systems - methods: n-body simulations 

1. Introduction 

Recent advances have brought the astronomical community to the verge of discovering 
a terrestrial planet in the habitable zone (HZ) of its parent star. Classically, the HZ is 
defined as the range of orbits for which a terrestrial mass planet (0.3 < M < 10 M®), with 
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favorable atmosp heric conditions, can sustain liquid water on its surface (IKasting et al.lll993 



Udrv et al. 


2007; 


Selsis et al. 


2007; 



20071 . for example), but no serious candidates for habitability have been found. 



Recently iBouchy et al.l (120081 ) used the HARPS spectrograph (Mayor et al 2003) to 
detect two planets around the nearby star HD 47186: a hot Neptune (HD 47186 b; 22.78 M^) 
with an orbital period of 4.08 days and a Saturn-mass planet (HD 47186 c; 0.35 Mjup) with 
a period of 3.7 years. The orbital elements for this system are given in Table [TJ A curious 
characteristic of the host star HD 47186 is that its mass and luminosity are ve ry similar to 
our own Sun, albeit with higher metallicity (Table 1. of Bouchy et al. ( 20081 )). Therefore 
the HZ of this system may be simila r to ou r own solar system's, and hence lies between 
the two known planets. iBouchy et al.l (120081 ) found that the combined radial velocity (RV) 
amplitude for this system is 4.30 m s~^ with 0.91 m s~^ residuals and reduced = 2.25. So 
potentially habitable planets may lie at or below the current RV detection threshold but may 
be discovered by further observations. Here we perform A^-body simulations to assess the 
stability of additional hypothetical planets within the HZ of HD 47186. The relatively small 
masses of the two known planets suggest an additional planet would likely also have low 
Masses less than 10 M^ are often considered the upper limit for terrestrial-like planets 



mass. 



dPoUack et a"l~lll996l : iBodenheimer et al. l[200ol ; iHubickvi et"al~ll2005h . lence an additional 



planet in the HZ of this system might have a mass conducive to habitability. 

As of Feb 27, 2009, 290 extra- solar planetary systems are known including 37 multiplanet 
system^. iBarnes fc Quinn (2004) found that many multiplanet systems tend to be close to 
the edge of stability such that slight perturbations to the system would result in destabi- 
lization jBarnes fc Quinn 2001; Gozdziewski et al. 2001; Gozdziewski &: Maciejewski 2001 



Erdi et al. II2004I ). This observation led to the "Packed Planetary systems (PPS)" hypothesis 
(Barnes & Raymond 2004; Raymond & Barnes 2005; Raymond et al. 2006; Barnes & Green- 
berg 2007; see also Laskar 1996) which postulates that if a region of stability is available in 
between two known planets, then that region will harbor an unseen planet (in other words, 
planetary systems tend to be "packed", with no large gaps between planets). To constrain 
the location of additional planets, a large number of studies have mapped out stable zones 
in known extra-solar planetary systems, mainly using mass less test particles (e.g., Rivera 
& Lissauer 2000; Jones et al 2001; Menou fc Ta bachnik 2003; Dvorak et al 2003; Barnes 
& Raymond 2004). Recently, iBean et al.l (120081 ) reported the d iscovery of a Saturn-mass 



planet in the narrow, well-characterized stable zone of HD 74156 (IBarnes fc Raymond! 12004 
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Raymond fc Barnesll2005l ). This constituted the first prediction of the orbit and mass of an 
extra-solar planet (Barnes et al. 2008a)@. 

The PPS model was furth er developed by the identificatio n of an analytic stability 
boundary in two-planet systems (IBarnes &: Greenberg]|2006l . 120071 ) . The boundary lies where 
the ratio of two quantities, /3//3crit, approximately equals unity. Here (3 depends on the 
total angular momentum and energy o f the system and /5crit depends only on the masses 
( iMarchal fc Bozislll982l : lGladmanlll993l ). If /3//3crit >1 (and no resonances are present), then 
the system is stable. Moreover, Barnes et al. (2008 a) found that systera s tend to be packed 
if /9//5crit ^1-5 and not packed when /5//3crit ^ 2. iBarnes fc Greenberg] (120071 ) pointed out 



that the vast majority of two-planet systems are observed with /9//3crit < 1-5 and hence are 
packed. 

Our investigation into possible habitable planets in HD 47186 is further motivated by 
the large separation between the two known planets: the P/Pcrit value for this system is 
6.13, the largest value among known, non-controversial systems that have not been affected 
by tidesif] This value is much higher than the packing limit, giving further support to 
the argument that at least one additional companion might exist between the two known 
planets. Moreover, the success of A^-body investigations in predicting the mass and orbit 
of additional planetary companions suggests a similar analysis on HD 47186 could yield 
similar results. But this time, the yet-to-be-detected planet could be habitable. In [|2] we 
describe our numerical methods and map ou t stable regions for a hy pothetical 10 Mq planet 
d between the two known planets. Following [Raymond et al.l (120081 ). we then investigate the 
perturbations of this hypothetical planet on the orbits of the known planets b and c. This 
approach allows us to exclude a portion of parameter space that is dynamically stable but 
offers a low probability of detection. Finally, in [|3] we discuss our results. 



2. Dynamical stable zones in HD 47186 
2.1. Numerical Simulations 



We started our simulations by assi gning the two know n planets from HD 47186 their 
best-fit orbits and masses from Table 1 (IBouchy et al.ll2008l ). and assuming mutual inclina- 



^It should be noted that HP 7 4156 d remains a controversial detection; Wittenmver et al. ( 2009h ii 
reanalysis of the lBean et al.l (j2008f ) spectra see no evidence for a third planet orbiting HD 74156. 



■^See http://www.astro.washington.edu/users/rory/research/xsp/dynamics/| for an up to date hst of 
/3//3crit values for the known extra-solar multiple planet systems. 
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tions of 1 degree. We also included a hypothetical planet d in between the orbits of planets 
b and c. From simulation to simulation, we varied planet ci's semi-major axis in the range 
0.06 < ctd < 2.32 AU in steps of 0.02 AU and eccentricity in the range < Cd < 0.6 in 
steps of 0.05, for a total of 1483 simulations. We assigned planet d a mass of 10 M^, which 
would induce in the star a radial velocity amplitude of 3.5-0.6 m s~^ in t he range that we 
considered (0.06 — 2.3 AU). We evolved each system with the Mercury ( IChambersI Il999l ) 
integrator, using the symplectic Wisdom- Holman mapping (Wisdom & Holman 1991). Each 
system was integrated with a time step of 1.0 day and conserved energy to better than 1 part 
in 10^. Simulations were stopped after 10 Myr or if there was a close encounter between any 
two planets of less than 3 Hill radii. Previous work has shown that most unstable planetary 
systems undergo instabilities on a timescale of 1 Myr or less (e.g., Barnes & Quinn 2004). 
Thus, our 10 Myr integrations should be sufficient to give a realistic depiction of the stable 
zone. 



2.2. Mapping stable zones 

In Fig. [1] we show the results of our simulations plotted in a-e space. The green and red 
filled circles represent stable and unstable simulations, respectively. The two known planets 
in HD 47186 are shown as blue filled circles, and the vertical bars represent the observational 
uncertainties in their eccentricities. The two curves overlaying the stable/unstable points 
are the crossing orbits for planets b and c. The blue shaded background is the "eccentric 
habitable zone" (EHZ) for this system (Barnes et al 2008b). The EHZ is defined as the 
HZ from Selsis et al (2007), assuming 50% cloud cover, and modified to take into account 
variations in the orbit-averaged fiux with eccentricity, which appears to be the key factor 
in determining surface temperature (Williams & Pollard 2002). Also labeled are some of 
the mean-motion resonances (MMRs) between the hypothetical planet d and outer massive 
planet c. Out of the 1483 simulations, 856 (57%) simulations were stable for 10 Myr, and 
the stable zone includes the inner half of the EHZ. Within this stable zone, several MMRs 
produce "islands" of instability associated with specific MMRs, notably 5d : Ic, 7d : 2c, and 
5d : 2c. Conversely, beyond the crossing orbit of planet c (black dashed curve overlapping 
EHZ), there are several MMRs, such as 2d : Ic, 7d : 3c (not shown, but lies in between 5d : 2c 
and 2d : Ic) and 3d : 2c that provide stability. This result is similar to that of Raymond et 
al (2008), who found regions of stability /instability associated with mean motion resonances 
in the 55 Cancri planetary system. 

Using just dynamical stability as a metric, all of the stable (green) points in Fig. [1] 
should be equally likely to support a 10 M® planet. To narrow down the possible orbit of 
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planet d, we cal c ulate a quantity called "Fraction of Time spent on Detected orbits" (FTD, 
Raymond et al. ( 20081 )). which takes into account the back-reaction of planet d on the orbits 



of planets h and c. The FTD provides a measure of the probability of finding planets h and 
c on their current best-fit orbits, including observational uncertainties (see Table 1), taking 
into account perturbations from planet d. In practice, the FTD simply evaluates the time 
spent by each observed planet {h and c) on their best-fit orbit, to within the uncertainties. 
A stable hypothetical planet can cause oscillations in the known planets' orbits that may be 
larger than the uncertainties. If the FTD value is small, then the oscillation amplitudes are 
large and the chance of observing the system in it's current configuration is small. Hence, it 
is less likely for a hypothetical planet d to have those orbital parameters. If the FTD value 
is close to 1, then it is more probable for planet d to be in that orbit. Note, however, that 
some systems may be observed in unlikely configurations (Veras & Ford 2008). 

Figure [2] shows the FTD for the stable region mapped out in Fig. [TJ yellow is the most 
likely region for planet d to exist (high FTD), and black is the least likely (low FTD; see 
color bar). Comparing Figs. [T]&[2], the region a < 0.2 AU is an unlikely location for planet 
d to exist. On the other hand, beyond a ~ 0.3 AU, the FTD value is uniformly close to 1. 
Therefore planet d does not significantly perturb planet c (otherwise there would have been 
low FTD values), and acts effectively as a massless particle. This leaves a large contiguous 
stable region with high FTD values from 0.3 — 1.3 AU. We see that planet d might exist in 
the habitable zone of HD 47186. 52.5% of the stable simulations are within the EHZ, but it 
is not clear how likely it is that planet d should exist there. 

We also performed 30 simulations to see if two additional planets could exist in the 
stable zone. We placed two planets of 10 each in the range 0.3 < ctd < 1-3 AU and 
< Cd < 0.6. We found that 16 cases resulted in unstable simulations and 14 stable. Most 
of the stable simulations (13, ~ 93%) had low eccentricity (e < 0.3). In some stable cases, 
both hypothetical planets were in the EHZ. 

Additional simulations assuming that the masses of the known planets b & c are twice 
their minimum mass, equivalent to an inclination angle of 30 degrees, did not show any 
significant change in the stability boundaries. We found that the inner stability boundary 
near planet b is not at all affected, whereas the region between 1.5-1.8 AU is slightly modified. 
Specifically, regions stabilized by 2:1 mean motion resonance near 1.5 AU, are now smaller 
with some of the stable points turning into unstable. Overall, there is no significant change 
in the stability of the HZ. 
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Discussion 



We have shown that a hypothetical 10 planet can exist between planets HD 47186 
b and c for orbital distances of roughly 0.1 - 1.3 AU (Fig. [1]). Our FTD analysis shows 
that this planet is more likely to be found between ~ 0.3 and 1.3 AU (Fig. [2]). We also 
find that two 10 planets with low eccentricities can exist between planets b and c. A 
10 planet between 0.3 and 0.7 AU (the high FTD region in Fig. [2]) produces an RV 
amplitude of 1 — 2 m/s, and should be detectable in the near future. It is possible that 
the two known planets in HP 47186 may have mig rated to their present location. Previous 
studies ( iRaymond et al.ll2006l : iMandell et al.l 120071 ) have shown that migrating giant planets 
don't preclude the existence of terrestrial planets. In fact, if the outer planet migrated 
inward then the most likely place for an extra planet would be just interior to the 2:1 MMR, 
which lies in the middle of the EHZ. Planets of ~ 10 or less are thought to be the largest 
potentially terrestrial planets with rocky surfaces and thin atmospheres. This provides an 
effective upper limit for habitability as we know it. Therefore HD 47186 offers the chance to 
detect a terrestrial planet in the HZ of a nearby, solar-type sta r. Roughly 30% o f systems 
are thought to harbor super-Earth sized planets (Mp < 10 M®, (IMayor et al.ll2009l )). If this 
frequency increases for lower-mass planets, the (non-)detection of HD 47186 d could shed 
light on the fraction of stars with terrestrial-like planets in the habitable zone, sometimes 
called 7]^. If, in a complex system such as this one, a terrestrial-scale planet can form and 
survive, then perhaps i]^ is relatively large, increasing the likelihood the spacraft such as 
Kepler and Terrestrial Planet Finder will find habitable planets. We encourage observers 
to use our results to assist in the search for HD 47186 d. 



The HD 47186 system is also another critical test of PPS hypothesis, which asserts 
that planetary systems should not contain large empty gaps. The PP S hypothesis follows 
naturally from the planet-planet scattering model (IRaymond et al.ll2009l ) and clearly predicts 
that there exists an additional planet in the gap between planets b and c. In fact, two 10 
planets could exist in the stable zone if they are at low eccentricity. The connection between 
the simple, easily calculated value /9//3crit niight provide some insight into the number of 
planets that can exist between two known planets. It would be interesting, then, to see if 
this one parameter can be used to test or not whether a system is packed. 
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Fig. 1. — The stable zone for a hypothetical 10 M® in between the two known planets of 
HD 47186. The red and green circles represent stable and unstable simulations, respectively. 
The known planets HD 47186 h and c are indicated by blue filled circles: the error bars 
represent the observational uncertainties in the best-fit eccentricities. The curves on top of 
the stable/unstable zones represent s crossing orbits for p lanets h and c. The blue background 
shows the eccentric habitable zone ( Barnes et al.ll2008bl ) for this system. Also shown are the 
several MMRs between planet d and c. 
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Fig. 2. — The FTD ("Fraction of Time spent on Detected orbits") ( (IRaymond et al.ll2008l )) 
for the stable zone of HD 47186. The FTD value is shown in the color bar: yellow indicates a 
high likelihood of detecting planet d with those orbital parameters. The long gaps correspond 
to unstable zones arising from MMRs. 
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Table 1. Orbital parameters of planets around HD 47186 fiBouchy et al.ll2008l ) 



Parameters HD 47186b HD 47186c 

P [days] 4.0845 ± 0.0002 1353.6 ± 57.1 

e 0.038 ± 0.020 0.249 ± 0.073 

msin i (Mj^p) 0.07167 0.35061 

a (AU) 0.050 2.395 



